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Abstract of EP07201 25 

The system formes an image (SG) from a source 
image (M) using a shooting system (1), data 
capture unit (2) and numerising unit (3). The 
image becomes distorted by the shooting system. 
A numerical system of image processing (4, 47) 
conrects the distorted image (SG) to a corrected 
target image (TG). The Image processing system 
(4, 47) has a first sub-assembly (4) which 
determines the optical centre (OCn) of distortion 
from a test source and also a polynomial function 
of correction (FNopt/NLast). A second sub- 
assembly (47) has a memory to stock these data 
and a calculation block to apply them. 
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(54) Dispositif de formation dimage et procede pour effectuer des corrections de distorsions 
optiques geometriques dans une image 

(57) Le dispositif connprend des systemes de prise 
de vue (1), de saisie (2) et de num^risation (3) d'une 
image source (SI) distordue, un systems de traitement 
d'image numerique (4,47) de construction d'une Image 
cible (Tl) corrlgee en distorsions, correspondent a ('ima- 
ge source (SI), lequel inclut un premier sous-ensembie 
(4) de determination prealable ^ partir des donnees 
d'une image source TEST distordue (SG^), d*un centre 
optique (OC") et d'une fonction polynomiale (F") de cor- 
rection des distorsions radiales autour de ce centre op- 
tique et un second sous-ensemble (47) qui applique la- 
dite fonction polynomiale de correction ^ chaque adres- 
se de pixel de I'lmage cible (Tl) pour fournir I'adresse 
d'un point dans I'image source distordue (SI) ou trouver 
une donnee d'intensite ^ appliquer^ i'adresse de depart 
dans I'image cible. 

Le procede comprend le calcul, dans une image 
source TEST {SG^), d'un centre optique (OC") au centre 
des distorsions optiques, et d'une fonction polynomiale 
(F") de correction des distorsions pour construire une 
image cible TEST (Tl). 

Application : systemes d'imagerie vid6o ou numeri- 
que 
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Description 

L'inventlon concerne un dispositif de formation d'image comprenant : 

s un systeme de prise de vue, 

un systeme de saisie des donnees d'une image source distordue par le systeme de prise de vues, 

un systeme de numerisation de ladite image source incluant une premiere memoire d'image pour le stockage des 

donnees d'intensite de chaque pixel repere par une adresse dans une matrlce bidimensionnelle, 

un systeme de traitement d'image numerique construction d'une image cible corrigee en distorslons, correspondant 

10 a I'image source. 

L'invention concerne egalement un procede pour effectuer des corrections de distorsions optiques geom§triques 
produites par un dispositif de prise de vue dans une image. 

L'invention trouve son application dans la correction des images num^risees dans le domaine des systemes a R- 
^5 X, ou bien des systemes video, ou bien des systemes d'imagerie numeriques. 

Des distorsions geometriques significatives peuvent etre introduites par I'objectif d'une camera dans I'image pro- 
duite par cette camera, tout partlcuH6rement, si cet objectif est du type grand-angle. Ces distorsions geometriques 
sont le plus souvent du type dit en tonneau, ou du type dit en coussin. Ces distorsions geometriques apparalssent 
mgme si I'objectif de la camera est de tr^s bonne quality. 
20 if est deja connu de I'etat de la technique par la demande de brevet GB 2 256 989 un dispositif pour r§aliser une 

compensation des imperfections optiques dans une camera de television. 

Ce dispositif comprend une camera de prise de vue d'enregistrement de I'image formee par le systeme optique, 
et des moyens pour numeriser cette image qui incluent le stockage des donnees d'intensite de chaque pixel courant. 
Ce dispositif comprend en outre un correcteur d'erreur pour compenser les erreurs geometriques, d'enregistre- 
25 ment, et chromatiques et compenser ainsi les imperfections optiques du systeme de prise de vue a lentille. 

Ce correcteur d'erreur est command e par un control eur de correction. Ce controleur de correction regoit, via une 
interface, des informations exterieures, appropriees a le programmer pour qu'il foumisse au correcteur d'erreur des 
signaux de controle qui prennent en compte des parametres du systeme de prise de vue. Dans ces conditions, le 
correcteur d'erreur est capable de corriger les donnees des pixels. 
30 Dans le correcteur d'erreur, une memoire cartographique regoit du controleur de correction ces signaux de controle 

dependant des parametres du systeme de prise de vue. Cette memoire cartographique est tabulee en fonction de ces 
parametres de maniere a definir la cartographie qu'il est necessaire d'utiliser pour corriger les defauts du systeme 
optique avec un type de lentille donne et dans des conditions de prise de vue donnees. 

La sortie de la memoire cartographique est fournie a un interpolateur utilise pour ameliorer la definition de I'image 
35 de sortie. 

La memoire cartographique est calibree en fournissant a la camera de prise de vue une MIRE de test ayant un 
dessin de grille carree reguliere, et en ajustant manuellement les parametres stock6s, de maniere ^ obtenir en sortie 
du dispositif une grille corrigee. Cette operation peut etre realisee en affichant par exemple la sortie du dispositif sur 
un ecran muni d'un graticule superpose. 

40 Le dispositif connu du brevet GB 2 256 989 n'opere done pas les corrections des imperfections du systeme de 

prise de vue de maniere automatique. II faut lui fournir des donnees sur la distance focale de la lentille. la distance de 
prise de vue, et le taux de zoom utilise. II faut en outre entrer dans une memoire cartographique, pour chaque type de 
lentille et chaque condition de prise de vue, des tableaux dans lesquels seront puisees les infomnatlons relatives aux 
nouvelles adresses a affecter aux pixels de I'image cible, pour remplacer les adresses d'origine de ces pixels dans 

45 I'image source c'est-a-dire dans I'image directement issue du systeme de prise de vue et entach6e d'imperfections. 
Ainsi I'image cible peut etre corrigee des distorsions de I'image source. 

Aucun moyens automatiques de tabulation de la m6moire cartographique, ou moyens automatiques pour calculer 
des fonctions de correction a appliquer aux pixels de I'image source, pour obtenir les pixels de I'image cible ne sont 
enseignes par le document cite. Seuls des moyens de calibration manuels de la memoire cartographique sont dlvul- 

50 gues, incluant la superposition de Timage cible distordue d'une grille sur un graticule de reference et la correction par 
Taction manuelle sur des parametres d'entree du systeme de prise de vue. 

Un but de la pr6sente invention est de proposer un dispositif, ainsi qu'un procede, pour calculer des fonctions de 
correction des distorsions optiques d'un systfeme de prise de vue, pour obtenir des donnees de pixels permettant de 
construire une image cible corrigee. 

55 Un but de la presente invention est de proposer un tel dispositif et un tel proc6d6 pour calculer automatiquement 

de telles fonctions de correction sans avoir ^ tenir compte des parametres du systeme de prise de vue. 

Ces buts sont atteints au moyen d'un dispositif tel que defini dans le pr6ambule de la partie introductive, dans 
lequel, le systeme de traitement d'image inclut : 
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un premier sous-ensemble, de determination pr^alable k partir des donn^es d'Image d'une image source TEST 
distordue, d'un centre optique et d'une fonction polynomiale de correction des distorsions radiales autour de ce 
centre optique, 

un second sous-ensemble incluant une merrioire ;de*stockage du centre optique et de la fonction polynomiale 
5 determines par le premier sous-ensemble, un bloc de calcul qui applique la fonction polynomiale de correction k 

chaque adresse de pixel de I'image cible pour fournir I'adresse d'un point dans Timage source distordue ou trouver 
une donnee d'intensite a appliquer ^ I'adresse de depart dans I'image cible, et une seconde m§moire dMmage pour 
le stockage des donnees de I'image cible reconstruite. 

10 Un precede pour effectuer des corrections de distorsions optiques produites par un systeme de prise de vue, dans 

une image, comprend des etapes pr^liminaires de : 

a) saisie des donnees et numerisation d'une Image source TEST incluant le stockage des donnees d'intensite de 
chaque pixel, 

IS b) estimation d'un centre optique localise au mieux au centre des distorsions optiques de I'image source TEST et 

report de ce centre optique calcule dans une image numerique a construire, referencee Image TEST cible, repre- 
sentant I'image source TEST corrigee des distorsions optiques, 

c) estimation d'une fonction polynomiale pour faire correspondre a I'adresse d'un pixel dans I'image TEST cible 
corrig6e, I'adresse d'un point correspondant note point de reference, dans I'image source distordue, en se fondant 
20 sur I'hypothese de base que les distorsions optiques geometriques dans I'image source sont radiales autour du 

centre de distorsron, de sorte que lesdits points correspondants sont alignes au mieux avec le centre optique 
estime. 



2S 



35 



so 



Un tel procede peut en outre comprendre dans les etapes preiiminaires, des etapes iteratives de : 



d) estimation de la meilleure fonction polynomiale F^opt capable de minimiser, a une iteration d'ordre Nopt, I'erreur 
de modelisation realisee sur la localisation des pixels de I'image test cible construite a une iteration precedente 
selon I'hypothese des distorsions radiales autour d'un centre optique OC^ estime k ladite iteration precedente, 

e) estimation, h une iteration d'ordre NIast. d'un centre optique modifie optimise OC^iast et reestimation dans ces 
30 conditions d'une nouvelle fonction polynomiale F^^PKNiast) capable de minimiser davantage I'erreur de modelisation 

realisee sur la localisation des pixels de I'image cible reconstruite selon I'hypothese des distorsions radiales autour 
dudit meilleur centre optique OC^^^^ estime. 



L'invention est decrite ci-aprfes en detail, en reference avec les figures schematiques annexees dont : 



la FIG.1 A qui represente un dispositif de formation d'images incluant un systeme de prise de vue et deux sous- 
ensembles d'un systeme de traitement d'images numehques. Tun pour determiner et I'autre pour appliquer des 
lois de corrections geometriques des distorsions optiques du systeme de prise de vue ; 

les FIG.1 B et 1 C qui illustrent I'etape de calibrage pour I'acquisition d'une mire de test par le systeme de prise de 
40 vue ; 

la FIG. 2 qui represente un diagramme schematique de blocs formant le premier sous-ensemble du systeme de 
traitement d'image, pour mettre en oeuvre un procede de determination de lois de corrections etde determination 
du centre de distorsion, pour corriger geometriquement les distorsions optiques provoquees par le systeme de 
prise de vue ; 

^5 - les FIG.3A a 3G qui representent des images source et cible dans les differentes etapes du procede ; 

la FIG.4A qui represente une intersection dans la grille source et illustre le filtrage non lineaire pour extraire les 
points de reference ; 

la FIG.4B qui represente une intersection dans la grille source pour la determination du point de reference Fig le 
plus proche du centre de distorsion, par un critere d'allgnement ; 

la FIG.4C qui represente un sous-ensemble horizontal HSj de points de reference pour la determination du pas 
de la premiere grille TG^, par un premier filtrage median ; 

la FIG.4D qui represente un sous-ensemble central S^° de points de reference pour la determination du pas de la 
premiere grille theorique TG^, par un second filtrage median ; 

la FIG.4E qui illustre la formation de couronnes rectangulaires de taille croissance, pour la formation des couples 



pO 



les FIG.5A et 58 qui representent des images sources distordues en coussin et en tonneau respect ivement ; 
ia FIG.6A qui illustre la mise en correspondance d'un point de grille TR^^^^^^ et d'un point de reference dans 
une couronne rectangulaire a, et la FIG. 68 qui illustre la meme operation dans une couronne adjacente plus 
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grande p ; 

la FIG.7A qui illustre la determination des erreurs A[J d'alignement des couples PjJ sur le centre optique OC". ces 
erreurs etant dues a la distorsion ; 

la FIG.7B qui illustre la detennination des couples de rayoiis XJJ et Y[J pour {'estimation de la fonction polynomiale 

5 F" ; 

la FIG.7C represente les points des couples p^^^ idealement alignes avec le centre OC^^^^* qui minimise au 
mieux I'erreur de modelisation ; 

la FIG.8A qui illustre la determination d'un polynome F" pour relier au mieux les points formes par les couples de 
rayons et YJJ report^s respectlvement en abscisse et en ordonnee sur un graphe dans un repere gradu6 en 
10 pixels d'origine OC" ; 

la FIG.8B qui illustre la determination du premier polynome F^ a I'iteration zero dans la methode iterative et la FIG. 
8C illustre la determination d'un polynome ameliore F"^ a une iteration n ulterieure ; 

la FIG.8D qui montre les points des couples p^'^^* idealement places sur la courbe representative du meilleur 
polynome f'^'^s*<'^°p*) pour minlmiser I'erreur de modelisation. 

IS 

I - DtSPOSITIF (Fig.1 A. Fig.2) 

D'une maniere generale, un systeme de prise de vue de camera peut Inclure un systeme optique a lentille de type 
grand angle, ou bien un systeme de lentilles de type zoom qui produisent des distorsions optiques geom6triques dans 

20 I'image formee par cette camera. Ces distorsions apparaissent meme si le systeme optique est de tres bonne qualite. 

En reference avec la FIG.1 A. un dispositif de formation d'image comprend un systeme de prise de vue 1 ayant un 
systeme optique k lentilles, qui fournit une image distordue d'une scene. Ce systeme optique 1 a lentille est attache 
a une camera 2 qui convertit I'image optique en signaux electriques grace a un element a CCD par exemple. Ces 
signaux electriques sont fournis a un systeme de numerisation 3 capable de stocker les donn^es d'image numerisees 

2S c'est-a-dire les donnees d'intensite relatives a chaque adresse des pixels de I'image reperes dans une matrice bidi- 
mensionnelle, dans une premiere memoire d'images. Ces donnees d'image numerisee sont relatives a I'image distor- 
due. 

Ce dispositif de formation d'images comprend aussi un systeme de traitement de signal numerique qui comprend 
un sous-ensemble 47 qui traite les donnees d'image numerisees stockees dans la memoire d'image du systeme de 
30 numerisation 3 pour fournir des donnees d'une image numerique correspondante, reconstruite avec correction des 
distorsions, ces donnees etant finalement stockees dans une seconde memoire d'image. 

On appellera ci-apres image source SI les donnees d'image numeriques en sortie du systeme de numerisation 3 
relatives a une image distordue par le systeme de prise de vue, et image cible Tl les donnees d'images numeriques 
fournies par le sous-ensemble 47 du systeme de traitement d'image, relatives a I'image reconstruite corrigee en dis- 
35 torsion. 

Le sous-ensemble 47 du systeme de traitement d'image comprend. stockees dans une puce ou dans une table 
(LOOK UP TABLE) des lois de correction des distorsions pour construire I'image cible corrigee Tl. Ces lois dolvent 
etre determinees au moyen d'un procede comprenant des etapes preliminaires. Un des buts de invention est de 
fournir un tel procede pour determiner ces lois de corrections de maniere automatique. Les lois de corrections sont 
40 determinees une lois pour toutes dans un premier sous-ensemble 4 du systeme de traitement. Ensuite^ elles sont 
appliquees automatiquement au moyen du sous-ensemble 47 du systeme de traitement d'image appele second sous- 
ensemble 47. L'interrupteur INT sur la FIG.1 A permet de mettre hors circuit le premier sous-ensemble 4 du systeme 
de traitement, lorsque ces lois ont ete determinees. 

On decrit ulterieurement le procede de determination des lois de correction mis en oeuvre dans le premier sous- 
es ensemble 4 du systeme de traitement d'image. Ce premier sous-ensemble 4 pourra aussi §tre une extension du sys- 
teme de numerisation 3, ou etre integre dans ce systeme de numerisation 3. 

Dans le dispositif de formation d'image selon I'invention, tel que repr6sente sur les FIG.1 A et 2 : 

a/ Dans le premier sous-ensemble 4 du systeme de traitement, on met en oeuvre : des etapes de calcul pour 
so determiner une fois pour toutes : 

* des lois de correction geom^triques des distorsions produites par le systeme 1 de prise de vue, 

* et la localisation precise du centre des distorsions geometrtques, qui est aussi le centre optique du systeme 
a lentilles de la camera, dans la plupart des cas. 

ss 

b/ Dans le second sous-ensemble 47 du systeme de traitement, on met en oeuvre des etapes de calcul, pour 
construire. aux moyens des lois precedentes et de la connaissance du centre de distorsion, pour chaque image 
source SI distordue, une image cible Tl corrigee en distorsion. 
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* une methode pour determiner ces lois de correction jdes,.dJstors ions geometriques et pour determiner de iagon 
concomitante ta localisation precise du centre de distorslon, mise en oeuvre dans le premier sous-ensemble 4. 

5 • une methode pour construire I'lmage cible Tl en utillsant ces connalssances, mise en oeuvre dans le second sous- 
ensemble 47. 

11/ METHODE POUR DETERMINER LES LOIS DE CORRECTION ET LE CENTRE DE DISTORSION (RG.1, bloc 
4 : FIG.2, blocs 41 a 46) 

10 

Un procede pour determiner de maniere concomitante ces lois de correction et les coordonnees precises du centre 
de distorsion est decrit ci-aprds, en reference avec la FIG.2 qui en represents les differentes etapes par un diagramme 
schematique de blocs. 

Ce procede ne necessite aucune mesures preiiminaires telle que la mesure de la distance locale, ou la mesure 
IS de la distance de prise de vue, ou la mesure reelle du pas d'une grille test. 

Ce procede repose sur I'utllisation de Thypothese que tes distorsions geometriques engendrees par le systeme 
de prise de vue sont radiales, c'est-a-dire que, dans I'lmage cible construite par le second sous-ensemble 47 du sys- 
teme de traitement, un point corrige doit se trouver sur un rayon determine par le point distordu et par le centre de 
distorsion. C'est pourquoi ce proc6d6 necessite, en plus de la determination des lois de correction, la determination 
precise concomitante du centre de distorsion. Neanmoins ce procede ne necessite aucune mesure reelle precise 
preiiminaire, Dans ce procede, on ne cherchera pas a mesurer precisement, par un senseur par exemple, la localisation 
du centre optique de la camera, puisque seul le centre de distorsion de I'lmage source est interessant, du fait que 
les distorsions sont radiales autour du centre de distorsion. 

En reference avec la FIG.2, ce precede mis en oeuvre dans le premier sous-ensemble 4 comprend au molns les 
2S etapes sulvantes : 

II. A/ ACQUISITION D'UNE MIRE 

Cette etape comprend les sous-etapes suivantes. 

30 

II.A1/ Construction d'une mire (FIG.1 A) 

On realise une mire de test M. 
- On choisit de preference k cet effet de realiser le dessin d'une grille sur un support plan rigide. Les mailles de 
35 cette grille peuvent etre carrees ou rectangulalres. Dans un exemple, cette mire peut etre un dessin sur un support a 
fond blanc d'1 m x 1 .50m, representant des barreaux noirs horizontaux et verticauxformant des mailles carrees comme 
montre en M sur ta FIG. 1 A. 

II.A2/ Calibragede la camera (FIG.1B, FIG.1C, FIG.3A) 

40 

On place la camera face au plan de la mire M pour faire I'acquisition d'une image nette de cette mire. Un calibrage 
est realise avant d'effectuer la prise de vue. Ce calibrage consiste k : 

* rendre le support plan de la mire M perpendiculaire ^ I'axe optique X'X de la camera de prise de vue, comme 
45 illustre par la FIG. IB ; 

* rendre les barreaux de la grille paralleles aux lignes et aux colonnes de i'element a CCD de la camera 2, comme 
Illustre par la FIG.1C ; 

* placer le plan support de la grille de telle sorte que tout le plan Image de la camera solt couvert par des mailles 
de grille, comme illustre par ta FIG.3A. 

so 

Les conditions de calibrage consistent done essentiellement a verifier trois conditions d'orthogonalite, deux de ces 
conditions d'orthogonalite etant verifiees lorsque I'axe optique X'X de la camera est perpendiculaire au plan de la mire, 
et la troisieme lorsque les barreaux de la grille de mire sont paralleles respectivement aux lignes et colonnes ortlno- 
gonales de I'element k CCD de la camera. 
55 L'operatlon de calibrage peut etre effectuee en utillsant simplement des equerres. Ce calibrage n'est pas contrai- 

gnant car ii apparait que le procede selon I'invention est tres robuste a d'eventuelles imperfections de ce calibrage. 
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II.A3/ Acquisition d'une image de la grille servant de mire (FIG.1 A et FiG.3A), referencee grille source 



A raide de la camera 2, munie de son systeme optique a lentilles 1 , on effectue une mise au point et une prise de 
vue de la grille servant de mire M. Le systeme de numerisatlon 3 en fournrt une image num^risee sous forme d'une 
s grille, appelee grille source SG"^ montrant des distorsions optique comme representee sur la FIG.3A. 

Les distorsions peuvent etre en forme de coussin comme sur la FIG.5A, ou en forme de tonneaux comme sur les 
FIG.5B et 3A. 

Il.b/ EXTRACTION DE POINTS DE REFERENCE (FIG.2, bloc 41) 

10 

Une image de la grille source distordue est numerisee et stockee dans les moyens de numerisation et slockage 
3 montres sur la FIG. 1 . 

En reference avec la FIG.3A, dans I'image de la grille source SG'^ numerisee, on appelle ci-apres "points de 
reference" R|^ des points d'intersection des barreaux verticaux et horlzontaux de cette grille source distordue SG'^. 
IS Par horizontal et vertical, on entend paralleles aux lignes et aux colonnes de I'element a CCD de la camera, comme 
II resulte du calibrage. 

En ce qui concerne les notations : 

* pour les points de reference R^, k est un indice qui permet de num^roter les points de reference de la grille source ; 
20 ces numeros peuvent etre 1 a 500 par exempte, 

* pour les images de la grille source, par exemple SG'^, I'exposant permet de numeroter I'etape du traitement d'ima- 
ge. Ici I'exposant (-2) indique que cette etape est la deuxieme avant celle qui foumira I'image de grille source SG® 
des lois et du centre recherches. 

2^ En reference avec la FIG.2, la chaTne des operations formant les etapes du precede seton Tinvention, comprend 

une premiere etape mise en oeuvre dans le bloc 41 , consistant en {'extraction des points de reference l\ de I'image 
de la grille source SG"2, par la determination de leur adresse dans cette image source au moyen des sous-etapes 
suivantes : 

30 * un filtrage de I'image source pour rehausser {'intensity des zones d'intersection, 

* un seuillage de I'image filtree ; 

* une labellisation des points seuill6s pour former des zones, et une extraction du barycentre de chaque zone. 
ii.B1/ Rehaussement du niveau d'intensite des zones d'intersectlon par filtrage (FIG.3B : RG.2, bloc 41 > 

35 

Dans le cas pris comme exemple ou les barreaux de la grille mire sont representes en noir sur fond blanc ou clair, 
les zones d' intersection des barreaux verticaux et des barreaux horlzontaux sont noires comme montre sur la FIG.3A. 
Dans cette sous-etape. on realise un traitement de I'image numerisee de la grille distordue SG-^ pour rendre ces 
zones Claires, et pour rendre plus sombre tout le reste de I'image. On obtient ensuite I'image d'une grille-source 
40 notee SG"^ comme montre ensuite sur la FIG.3B. 

Ce traitement peut etre mene a bien par une premiere methods connue de I'homme du metier consistant en un 
filtrage lineaire au moyen de filtres lin^aires correspondant a des masques de correlation ayant la forme d'une croix. 

On aurait pu, aussi, au lieu d'une mire de test en forme de grille, choisir une mire de points. L'experience a montre 
que I'enregistrement par la camera d'une mire de points est moins bon que celui d'une mire en forme de grille, a cause 
4S d'eftets d'eblouissement. L'enregistrement d'une mire en forme de grille a en outre I'avantage de fournir davantage 
d'informations exploitables pour mener a blen le procede selon invention. 

On propose ci-apr6s une m6thode de traitement qui met en oeuvre un filtrage non Iin6aire notablement plus efficace 
que le filtrage lineaire connu. Les avantages de ce filtrage non lineaire par rapport k un filtre lineaire conventionnel 
sont qu'il : 

so 

* ne produit pas de tache d'eblouissement, 

* est tres robuste ^ la distorsion de I'intersection, 

* il est facile ^ mettre en oeuvre. 

55 Dans I'etape d'extraction des points de reference qui est op^ree sur INmage source numerisee SG~^ de la grille, 

le rehaussement des intensites des intersections des barreaux de grille est essentiel. II est done important de reaiiser 
un tres bon filtrage. 

Une intersection d'un barreau de grille horizontal avec un barreau de grille vertical est representee schematique- 
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ment sur la FIG.4A. Les bords fJu barreau horizontal et du barreau vertical sont representes par des lignes distordues 
puisque Ton se trouve dans Timage source distordue. 

Dans cette etape, on passe un fillre en chaque point de I'image source SG'^ de la FIG.3A. Par ce filtrage. on 
determine : 

5 

* C le centre du filtre, 

* S, N, E, W quatre points cardinaux, 

* SE, SW, NE, NW quatre points en diagonale. 

10 Les 8 points entourant le centre C du filtre sont enti§rement d§finis par des distances d1 et d2. 

* d1 et la distance mesuree en pixels entre le centre C et les points cardinaux (d1 = C-S, C-N, C-E, C-W), 

* 62 est la distance mesuree en pixels entre le centre C et les points en diagonale (d2 = C-SE, C-SW, C-NE, C-NW). 

IS La distance d1 est choisie telle que les points cardinaux retenus sont sltu6s k I'interieur du dessin des barreaux 

de grille quand le centre C est situe dans la zone d' intersection. Ce resultat est obtenu en choisissant di de I'ordre de 
1 a 5 pixels ; d'une maniere generate de I'ordre de la moitie de I'^paisseur d'un barreau de grille mesuree en pixels 
dans rimage distordue. 

La distance d2 est choisie de telle sorte que les points diagonaux correspondants sont situ6s dans le fond de 
20 I'image, c'est-a-dire en dehors des regions de dessin de barreaux de grille. Get effet est obtenu en choisissant d2 
mesur6 en pixels, de Tordre de la moitie du pas de la grille distordue dans I'image numeris§e. Par exemple pour une 
grille ayant un pas de I'ordre de 50 pixels, d2 pourra etre de 20 pixels. 

En choisissant des distances dl et d2 raisonnables que Thomme du metier peut determiner par quelques essais 
de routine, sans qu'aucune precision soit requise, ce filtrage permet d'obtenir que, lorsque le centre du filtre est ^ 
25 I'interieur d'une zone d'intersection, alors les points cardinaux S, E, W sont a I'interieur des regions des dessins des 
barreaux de grille, et les points en diagonale sont effectivement dans le fond, et ceci meme dans les regions de I'image 
distordue ou Ton trouve les plus fortes distorsions. 

Dans le cas exemple de la grille formee de barreaux noirs sur fond blanc, avant filtrage, dans I'image SG-2, rin- 
tensite de chacun de ces cinq points, constitues par le centre et les quatre points cardinaux, est faible, alors que 
30 I'intensite des quatre points diagonaux est elevee. 

Par le filtrage non lineaire selon I'invention, on realise une mesure exprimee par la mesure du critere FILTsuivant : 

FILT = Min(NW, NE, SE. SW) -Max(C, N. S, E. W). 

Dans ce critere Min(NW, NE, SW, SE) signifie que Ton lvalue le minimum d'intensite relatif aux points en diagonale, 
35 et Max(C, N, S, E, W) signifie que Ton evalue le maximum d'intensite relatif aux points cardinaux incluant le centre. 

Lorsque le filtre est correctement centre sur une intersection, normalement I'intensite de chacun des points en 
diagonale est grande done le minimum d'intensite evalue sur les intensite de ces points est encore grand. D'autre part, 
normalement I'intensite de chacun des cinq points, centre et points cardinaux, est faible, done le maximum d'intensite 
evalue sur les intensit^s de ces points est encore faible. 1 1 en resulte que la difference entre ce Min et ce Max est grande. 
40 Le filtre non lineaire realise revaluation de ce critere FILT en chaque point de I'image distordue de la FIG,3A. Le 

resultat du calcul du critere FILT constitue la sortie du filtre non lineaire selon I'invention. 

Les regions d'intersection detectees sont celles ou la mesure de ce critere est !a plus grande. L'image obtenue 
resultant de ce filtrage est une nouvelle grille source SG"^ et montre, comme sur la FIG.3B, que les regions d'inter- 
section sont rehaussees en intensite, et que les autres regions de I'image source distordue sont assombries. 

45 

M.B2/ Seulllage applique aux intensites des points rehausses (FiG.2, bloc 41 : FIG.3C> 

Apres la sous-6lape de rehaussement de I'intensite des points des zones d'intersection, on realise une sous-^tape 
de seuillage sur I'intensite de ces points. 
50 On propose ci-apres une methode de seuillage performante pour I'application a la detection des points de refe- 

rence. Selon cette methode, un seuillage est mene a bien dans I'image rehaussee SG*^ en selectionnant un nombre 
reference Nb representant le nombre des points de plus forte intensite. 

L'homme du metier qui a acquis I'image de la grille-source distordue SG"^ telle que representee sur la FIG.3A, 
effectue : 

55 

* le comptage du nombre N1 des intersections dans I'image, 

* le comptage, compte tenu de la numerisation, du nombre N2 approximatif de pixels contenus dans une zone 
d'intersection : chaque zone d'intersection a quatre cot^s, et sa surface est donnee par le produit des epaisseurs 
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en pixels des barreaux de grille horizontaux et verticaux dans I'image distordue, comme montre par exemple sur 

la FIG.4A. 

* le calcul du nombre Nb recherche egal au produit du nombre d'intersections N1 dans I'image distordue par le 
norDbre de pixels approximatif N2 dans une zone d'intersection. 

5 

Ce nombre Nb constitue le seuil 

Nb = N1 xN2 

A I'aide de ce seuil, on extrait les pixels appartenant aux zones d'intersection. Cette extraction consiste a retenir un 
10 nombre egal au nombre calcule Nb, parmi les pixels ayant les intensites les plus elev6es apres I'operation de filtrage. 

dans ('image rehaussee SG"V 

A cat effet, on realise un histogramme en comptabilisant, pour chaque niveau d'intensit^, le nombre de pixels 

ayant cette intensite. En partant du niveau d'intensite le plus 6leve, I'operation de seuillage consiste a sommer le 

nombre de pixels par niveau d'intensite decroissant, et a arreter la sommation lorsque le nombre Nb seuil est atteint. 
IS Ainsi les pixels de plus haute intensite sont retenus et leur donnees d'adresse et d'intensite sont memorisees. 

A Tissue de cette operation, I'image realisee a partir de la grille test n'est plus formee que de zones brillantes 

contenant un certain nombre de pixels dont I'intensite est superieure ou egale a I'intensite correspondant au seuil Nb 

sur rhistogramme, le fond etant uniformement sombre et depourvu des details de barreaux que Ton voyait encore 

apr6s filtrage sur I'image SQ-^ de la FIG.3B. 



20 



II. B3/ Etiquetage des points seuilles et calcul du barycentre des zones d'intersection (FIG.3C : FIG. 2, bloc 41) 



Les pixels extraits par seuillage a I'etape precedents sont regroup^s par une methode d'Stiquetage. Les zones 
d'intersection contenant chacune quelques pixels seuilles regroupes sont identifiees, puis leur barycentre est calcule. 
^s Par le calcul du barycentre, on ramene chaque zone d'intersection a un seul point, defini par un couple de coor- 

donnees qui ne sont pas forcement entieres puisqu'elles sont le resultat d'un calcul de barycentre d'une zone de 
plusieurs pixels. 

En reference avec la FIG.3C, les barycentres ainsi extraits constituent des points ci-apres appeles points de 
reference R^, ou k est le nombre qui permet de numeroter chacun des points de reference dans I'image source SG^ 

30 ainsi obtenue. 

A chaque R^. correspond une adresse dans I'image source SG^ maintenant obtenue. donnee par les coordonnees 
du barycentre calcule. Ces coordonnees sont notees rx^ selon Taxe des x et ry^ selon I'axe des y d'un repere de 
coordonnees fixe Qx, Cly utilise dans cette image source SG^. Ce repere est orthogonal et ses axes coincident 
respectivement avec une ligne et une colonne de pixels de cette image. 
35 Sur la FIG.3C, les barycentres ou points de reference extraits sont representes par des points noirs sur fond 

blanc pour la commodite du dessin, alors que dans I'exemple decrit, ils sont en realite blanc sur fond noir 

li.C/ ESTIMATION D'UNE PREiVIIERE GRILLE THEORIQUE fFiG.2, BLOC 42) 

40 Le precede de Tinvention comprend maintenant des etapes pour construire une image cible d'une grille appelee 

grille theorique TG" qui represente la grille SG° de I'image source (FIG.3C) corrigee de ses distorsions geometriques. 

A Tissue de ces etapes, le proced6 selon invention fournit des fonctions de transformation, permettant de cons- 
truire de maniere tres precise cette grille theorique TG" corrigee, a partir des donnees d'image de la grille source SG^ ; 
ces donnees d'image sont constituees au depart par les coordonnees rx,^, ry^ des points de reference R,^ dans le repere 

45 Qx, Qy. 

Selon I'invention, la grille theorique corrigee TG" est construite en se fondant sur I'hypothese que I'image source 
est distordue de maniere radiale autour d'un centre de distorsion par rapport k I'image cible ^ construire. 

Pour arriver a une construction precise de cette grille theorique amelioree TG" telle que representee sur la FIG. 
3G, le precede comprend d'abord une etape mise en oeuvre dans le bloc 42 de la FIG. 2, comprenant elle-meme des 
so sous-etapes pour la construction d'une premiere grille theorique TG*^ : 

* une estimation d'un point appel6 centre de grille cible, note GC^, localise a une intersection d'un barreau hori- 
zontal et d'un barreau vertical de cette premiere grille theorique. II est repere parses coordonnees en pixels gc>^, 
et gcyO dans un repere orthogonal Qx, Qy de la premiere grille theorique, dont I'axe des x est parallele aux 

ss lignes de pixels et I'axe des y est parallfele aux colonnes de pixels, et dont I'origine correspond a I'origine O du 

repere fixe dans I'image source SG^, 

* une determination de la position d'un centre optique OC^ considere comme un centre de distorsion approximatif 
pour cette premiere grille theorique TG°. 
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* une estimation d'un pas GJtJfi pour cette premiere grille thdorique, defini par les distances en pixels gmx^, gmyO 
entre les intersections de barreaux de cette premiere grille cible, selon Taxe des x et selon I'axe des y, dans le 
repere fixe Qx, Qy, 

* une mise en correspondance de points TR^ He la pj-eiriierfe grille th^orlque avec des points de r6f6rence F\ de la 
5 grille source distordue SG^. La localisation de ces points TR^ situes k Tintersection des barreaux de la premiere 

grille th§orique, est definie k partir des estimations du centre GC^ et du pas GM^ realises pour cette premiere 
grille th^orique. 

Sur les FIG.3C et suivantes lllustrant le present expose, or? a represents les points de reference R^ Timage source 
10 par des petits ronds, et les points d'intersection TR^ de la grille theorique par des petites croix en noir sur fond blanc, 
quel que soit leur niveau d'intensite reel pour la simplification des dessins. 

II. CI/ Estimation d'un centre GC^ pour la premiere grille theorique 

IS Cette premiere grille theorique est entierement definie des que sent determines un centre GC^ et un pas GM^. 

Les barreaux de la grille theorique sont orthogonaux et regullerement espaces. II taut noter que ce centre de la 
grille theorique pourra se trouver a n'Importe quel endroit de I'image, et meme, evenluellement en dehors de I'image. 

Dans cette etape, I'estimation d'un centre GC° pour la premiere grille theorique, est eftectue autour du centre de 
distorsion de la grille source distordue, pour mener ^ une estimation facile du pas GIVP de cette premiere grille thtorique 

20 jqo 

En appliquant I'hypothese des distorsions radiales, il apparait que le centre reel de distorsion de la grille source 
SG® est situe dans une region de I'image source ou les distorsions sont minimales. 

L'estimation d'un centre GC^ pour la premiere grille theorique comprend done la succession d'operations 
suivantes : 

25 

* une estimation d'un point de reference de depart R^ dans I'image de la grille source SG^, proche du centre de 
distorsion, 

* te report de ce point de depart R^ comme centre de grille GC^ et comme centre optique approximatif OC^ dans 
I'image de la premiere grille theorique TG^. Le centre optique de la premiere grille theorique est different du centre 

30 optique de la camera 1 . 2 de la FIG.1 A. 

11. C1 a/ Estimation d'un point de reference proche du centre de distorsion dans I'lmaqe source SG^ (FIG.2, 
bloc 42 : FIG.4B^ 

35 D'une maniere generale toute fonction capable de fournir l'estimation d'un point de reference, note R^, ,le plus 

rapproche du centre de distorsion de la grille source SG^, peut etre utilisee. On presente ci-apres une methode simple 
et appropriee a mettre en oeuvre cette operation. 

En reference avec la FIG.4B, un point de reference R^ de la grille source tel que represents sur ta FIG.3C est 
proche du centre de distorsion reel, torsqu'il est localise dans une region de cette grille source SG^ ou les distorsions 

40 sont les plus faibles. Sur la FIG.4B, on a represente un point de reference note R^ et un ensemble de points de 
reference voisins entourant ce point de reference F^^. L'ensemble des points peut etre decompose en deux sous- 
ensembles de points de reference : 

* un sous-ensemble horizontal not6 HS^^ compose de points de reference situ6s sur fe meme barreau de grille 
45 horizontal dans {'image source, que le point de reference de depart R^^, 

* un sous-ensemble vertical note VS|^ compose de points de reference situes sur le meme barreau de grille vertical 
dans I'image source que le point de r6f6rence de depart R^. 

On cherche a selectionner, parmi tous les points de r6f6rence de la grille source SG° celui rJ qui est le plus proche 
50 du centre reel de distorsion. 

Pour la determination des points de reference R,^, I'homme du metier avait dej^, dans I'etape de filtrage 1.B1/ 
illustr^e par la FIG.4A, affectum une determination approximative du pas de grille source, note GM-^ sur la F1G.3A. Ce 
pas approximatif de grille source est defini pour deux distances gmx*^ et gmy^ le long des axes de coordonnees Qx 
et Qy de I'image source SG^. 

ss Selon invention, on propose maintenant un critere d'alignement par rapport a une bande horizontale et une bande 

verticale de largeurs egales aux composantes horizontale et verticale du pas approximatif GM"^ c'est-^-dire respecti- 
vement de largeur gmx-^ et gmy ^ ces bandes etant centrees sur le point de reference teste R^^. Dans ces bandes, 
on s6lectionne les points formant le sous-ensemble HS^ et le sous-ensemble VS^. respectivement. 
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De rhypothese radiale, ii ressort que : plus le point de reference R^^ teste est prcx:he du centre de distorsion, plus 
les points des sous-ensembles VS^ et HS^ sont alignes, et plus nombreux sont les points de ces sous-ensembles 
dans chacune des bandes respectives de largeur gnnx"^ et gmy^ comme montre sur la FIG.4B. 

En reference avec ie schema de la FIG.4B> on propose le crrtdre d'alignement sulvant : 

5 

AL(ka) = 2 [(0,5 gmx'' - | rx^ - rx^J)2+(0.5 gmy'^ - |ry^ - i^yKaD^] 

relation dans laquelle nc^ ry^ sont les coordonnees d'un point de reference l\ de Tun des sous-ensemble HS^ et 
70 VS^ ; et dans laqueile rx^a, ry^a sont les coordonnees du point de reference test§ pour former le point de reference 
le plus proche du centre de distorsion. 

Le critere propose AL(ka) est calcule sur chaque ensemble de points de reference de la grille source SG® 
(FIG.3C). II apparait que : 

75 * plus la mesure representee par ce critere AL(ka) est grande, et plus il existe des points appartenant a cet ensemble 
. Done, plus il existe des points de reference de la grille source SG° suffisamment bien alignes avec le point 
de reference teste, pour ne pas depasser, avec ce point R^^ teste, une distance geometrique egale h la moitie de 
la composante du pas approximatif gmx-^ horizontalement, et gmy^ verticalement, 

* ce critere est done une description naturellement bonne d'une region peu distordue. 

Dans ce calcul du critere d'alignement AL{ka), on peut affecter les points de reference de I'ensemble de poids 
d'autant moins eleves que ces points sont plus eloignes du point de reference R^^ test6. 

II. CI b/ Report du point de reference retenu R j^ dans la premiere grille theorique a construire TG^ 

25 

Les points de reference R|^ de ia grille source SG^ de la FIG.3C sont testes en apphquant le critere d'alignement 
AL(ka) propose, et un point de reference RJ pour lequel ta mesure du critere est la plus grande est retenu. Ce point 
rJ est reporte dans I'lmage de la premiere grille theorique TG^ a construire. II est alors reference centre GC<> de la 
premiere grille theorique TG^ et see coordonnees gcx^ et gcy^ sont memorisees dans le repere Hx, Qy. On note 
30 que le centre de grille dans la grille theorique est toujours a une intersection de barreaux. 

Un centre optique OC° de la premiere grille theorique est aussi choisi en coincidence avec ce centre GC^. Le 
centre optique de la grille theorique est ci-apres defini comme le meilleur centre des distorsions connu dans Tetape 
du proced6 en cours. 

3S I.C2/ Estimation d'un pas GM^ pour la premiere grille theorique 

La valeur approximative GM'^ du pas de grille utilis^e par I'homme du metier dans les phases du procede du debut 
jusqu'a la sous-etape precedente d'esttmation d'un centre de grille GC^ pour la premiere grille theorique TG°, doit 
maintenant etre amelioree. A cet effet, on propose ci-apres une methode de calcui du pas de la premiere grille theorique 
40 effectue en deux operations de filtrage : 

* un premier filtrage pour determiner une valeur grossiere d'un pas GM*^ de la premiere grille theorique, 

* un deuxieme filtrage pour determiner une valeur plus precise d'un pas GM*^ pour la construction de la premiere 
grille theorique TG^. 

45 

II.C2a/ Premier filtrage pour determiner une valeur grossiere d'un pas GI\/l°o pour la premiere grille theorique 
(FIG.4C> 

Dans la sous-etape precedente, on a realise Testimation d'un centre de grille GC^ pour la construction de la pre- 
50 miere grille theorique, comme le point de reference R^ le plus proche du centre de distorsion dans I'image source SG^. 

En reference avec la FIG.4B, a ce point de reference rJ etait associe un ensemble de points de reference Sj, 
decomposable en deux sous-ensembles de points de reference : 

* un sous-ensemble horizontal note HS^ compose des points de reference situes sur le meme barreau horizontal 
55 que le point de reference R^ dans I'image source SG^, delimite par une bande horizontale de hauteur gmy** qui 

est la composante verticale du pas de grille approximatif GM"^; 

* un sous-ensemble vertical note VS|^ compose des points de reference situes sur le meme barreau vertical que le 
point de reference R^ dans I'image source SG^. delimite par une bande verticale de largeur gmx~i qui est la 
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composante horizontale du pas de grille approximatif GM~^ 

En reference avec la FIG.4C, on realise le premier calcui grossier des deux composantes gmx^^ et gmy^^ d'un 
pas pour la premiere grille thtorique : 

* en calculant la composante gmx^^ selon I'axe des Qx du rep6re fixe comme le r^sultat du filtrage par un filtre 
median applique sur les distances HD^ entre les points de reference du sous-ensemble horizontal HS^, 

* en calculant la composante gmy^ selon I'axe Qy du rep6re fixe comme le resultat du filtrage par un filtre median 
applique sur les distances entre les points de reference du sous-ensemble vertical HS^. 



A cet effet, on ordonne les points de reference dans le sous-ensemble horizontal HS^. On fait de meme dans le 
sous-ensemble VS^ non represente ; puis on calcule les valeurs des distances entre ces points de reference, et on 
ordonne ces valeurs de la plus petite a la plus grande. Ensuite, on applique un filtre median sur les valeurs ordonnees 
dans les ensembles HSj et VSj. horizontal ement et verticalement. Le filtre median retient la valeur situee au milieu 
75 des valeurs ordonn6es horizontalement pour gmx^° ; et de meme verticalement pour gmy°o. ' 

Les filtres medians ont pour avantage de tournir des resultats robustes aux erreurs de localisation des points 
traites. En particulier le resultat reste correct, dans la bande horizontale ou verticale de points traites, meme lorsqu'un 
point est manquant parce qu'il n'a pas ete extrait, ou bien lorsqu'un point extrait existe par erreur Les filtres medians 
donnent des resultats mellleurs que les catculs simples de moyennes. 
20 Les mesures effectuees dans cette premiere phase de filtrage constituent un couple robuste gmx^^, gmy^° pour 

une valeur grossiere du pas dans la construction de la premiere grille theorique. Cette valeur grossiere n'a cependant 
pas la precision suffisante pourconstruire la premiere grille theorique. C'est pourquoi, on precede ensuite k un second 
filtrage pour fournir une valeur du pas de grille suffisamment precise pour debuter cette construction. 

2B II. C2 b/ Deuxieme filtrage pour determiner une valeur plus precise d'un pas GM^ pour la construction de la 
premiere grille theorique (FIG. 2, bloc 42 : FIG.4D) 

En reference avec la FIG.4D, on selectionne un sous-ensemble S'^ de 8 points de reference situes autourdu point 
^ de reference retenu comme le plus proche du centre de distorsion dans I'image source SG®. La recherche de 
30 chacun de ces voisins est taite en utilisant la valeur du pas approximatif GM^o determinee dans le premier calcul. 
Chacun de ces 8 voisins est identifie comme le point de reference correspondant le mieux a ce premier calcul.' 

De nouveau, on applique un filtre median qui opere maintenant sur les valeurs des distances entre les points de 
reference ainsi definis. On utilise dans ce deuxieme filtrage 6 intervalles horizontaux HS° et 6 inten/alles verticaux 
VS° ; on ordonne ces inten/alles par valeurs croissantes et on determine par filtrage median les nouvelles valeurs des 
35 composantes horizontale gmx^ et verticale gmy^ du pas de grille, notee GM^, situees au milieu de ces valeurs ordon- 
nees en horizontal et en vertical respectivement. qui sont plus precises que les valeurs determinees par le filtrage de 
I'etape precedente, du fait que le filtrage est uniquement realise dans une region de I'image source proche du centre 
de distorsion, et done presentant les plus faibles distorslons. 

On note que, dans chacun de ces deux filtrages, des points de r6f6rence utilises pour passer les filtres medians 
40 peuvent etre errones ou manquants. Mais du fait que ces premier et second filtrages sont bas6s sur la selection des 
valeurs median es des intervalles consideres, les resultats ne sont pas notablement affectes. 
II en resulte que t'on dispose k la fin de cette sous-etape en deux operations : 

* d'une mesure robuste, 

45 * el d'une mesure suffisamment precise d'un pas GM^= (gmx° gmy*^), pour la construction de la premiere grille 
theorique. Ces donnees sont memorisees. 

Dans la methode selon I'invention, la robustesse des donnees d'image de depart pour construire la premiere grille 
theorique, k savoir : 

so 

* le centre de grille GC^ determine par les coordonnees d'un point de reference proche du centre de distorsion, 

* le centre optique OC^ egalement determine par R^ 

* le pas de grille GM^, 

55 sont tres importantes, car I'invention vise a fournir une methode automatique de determination des fonctions de cor- 
rection. 
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li.C3/ Estimation de premiers couples P {j formes de points de reference Ri,^de la grille source SG^ et de points 
TR ^ de la premiere grille theorique TQQ (FIG.2, bloc 42) 

5 Ayant maintenant la connaissance. 

* des coordonnees gcx^, gcy^ du centre de grille GC^, 

* des composantes gnrix®, gmy*^ du pas de grille GM^ dans le repere fixe Qx, Qy, on construit la premiere grille 
theorique en plagant, dans ce repere fixe, les intersections de ses barreau notes TR° , ou m est un indice qui 

10 permet de numeroter chaque point d'Intersection dans la grille theorique.. 

Dans cette etape, on travaille sur I'image source SG^ des points de reference telle que repr6sentes sur la FIG. 
3C, et sur la premiere grille theorique TG^ comma definie plus haut. En partant de ces deux Images SG^ et TG°, 
on determine pour la premiere fois des couples notes P^. de points correspondants formes d'un point de r6f6rence 
15 et d'un point de grille theorique TR^ . 

On ecrit cette correspondance sous la forme : 

= I^f(k)' ^'^g(k) ] 

ou f et g sont des fonctions de correspondance entre les points de reference de grille source et les points d'inter- 
20 section recherches TR^ de la grille theorique. En effet, rien n'Indique que le point TR^ trouve dans la grille theorique 
TQO pour correspondre a Rj^ dans la grille source SG^ aura le meme numero, et done le meme indice que le point de 
reference R|^. 

Par exemple Redans SG** pourrait avoir le numero 100 (k = 100) alors que TR^ dans TG° pourralt avoir le numero 
200 (m = 200). Dans ce cas, on ecrirait que f(k) = g(k), c'est-a-dire que f(100) = g(200) lorsque les points R,^ et TR^ 

25 sont apparles, et on determineralt ainsi une relation de mise en correspondance des points du couple P^. Ce pourrait 
etre par exemple f(100) = g(200) = 120. 

Dans la region de faible distorsion proche du centre de grille GCO, {'operation de premiere mise en correspondance 
consiste a s6lectionner les couples Pj formes chacun d'un point de reference R^^y.^ et d'un point TR^^j^^ de la premiere 
grille theorique, le plus proche I'un de I'autre par un critere de proximite simple. 

^o Au fur et a mesure que Ton s'eloigne du centre, cette fagon d'operer ne peut plus etre utillsee, car les points des 

couples deviennent de plus en plus eloignes I'un de I'autre, et vers les bords des images source SG^ et cible TG^ 
consid§r6es, on se trouve dans la situation ou le point de grille theorique doit etre mis en correspondance avec un 
point de reference selectionne parmi plusieurs points de reference qui semblent tous a priori de bons candidats. 
Dans le but de resoudre ce probleme, on part done du centre de grille GC^, et on forme le premier couple dans 

35 cette premiere grille theorique constitue du centre de grille GC^ situe a une intersection de barreaux, et avec le point 
de reference R^ qui coincide avec GC°, puis on construit des seconds couples Pj autour de ce centre de grille GC^. 
Ensuite, en s'appuyant sur les connaissances acquises de ces premier et second couples P^, on construit de nouveaux 
couples de proche en proche, en progressant pas a pas du centre de grille GC^ vers les bords par une methods decrite 
cl-apres. 

40 Selon cette methode, pour construire pour la premiere fois des couples P^ en progressant pas a pas du centre 

GC° vers les bords des images, a chaque pas, on d^finit une couronne rectangulaire de taille croissants autour du 
centre GC^. 

AInsi, dans a la fois la grille source SG^ et la premiere grille theorique TG^, avec le point de reference R^ et le 
centre de grille GG^ en coincidence : 



45 



so 



la premiere couronne se compose du centre de grille GC^ lul-m§me, 

la seconde couronne rectangulaire se compose de ses 8 voisins directs de GG^: 

la troisidme couronne rectangulaire se compose de ses 16 seconds voisins de GC^. 



Sur la FIG.4E, les points TR^^^ de la grille theorique sont representes par des croix quand ils sont suffisamment 
eloignes des points de reference Hj^k), eux-memes representes par des ronds, pour en etre distingues. 

Plus la couronne traitee est 6loignee du centre de grille GC^, et plus les couples sont formes de points distincts 
et distants. 

En reference avec la FIG.6A, dans une couronne rectangulaire courante a, apres la mise en correspondance d'un 

55 point de reference Rf(ka) d'un point TR^^j^^^ de grille theorique, on trace un rayon passant par le centre de grille GC^ 

et par le point de grille TR^ teste dans cette couronne courante a. Et on mesure le vecteur d'erreur : 
g(K(x) 

^ ka- ' (ka) • "(fka) 
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Ce vecteur V^^ constitue I'erreur de mise en cprrespondance des points du couple P^^ dans cette couronne a, erreur 

due a ta distorsion. 

En reference avec la FIG.6B, dans une couronne p suivante plus eloignee du centre de grille GC^, un point de 
^ r6f6rence est recherchd dans I'image source SG^ pour correspondre ^ un nouveau point test6 TR^^j^p^ de la premiere 

grille thtorique TG°. A cet effet, on a nnemorise les vecteurs d'erreur V^^ relatifs ^ tous las couples P^^ de la couronne 
pr6c6dente a. On trace le rayon passant par le centre GC^ et par le point test6 TR^^j^pj. On recherche dans la couronne 
^ precedente a le point de grille th6orique OTR^^j^^^ le plus prfes de ce rayon (GC®, TR^^^^^^ que Ton vient de tracer. Ce 
point de grille "^'R^^i^^j correspondait a un point de reference Rf^^ot) P^^ '® vecteur d'erreur 9V^,^^ defini plus haut. Ce 
vecteur d'erreur V^^ est alors applique au nouveau point test6 TR^^j^p^ pour estlmer le point de reference correspon- 
dant. Le point de reference Rf(icp) le plus proche de cette estimation est retenu s*il est trouve dans une fourchette de 
distance acceptable d6fini par les composantes du pas de grille th6orique gmx^, gnny^. Un nouveau couple P^^ est 

alors forme dans la mesure ou un point de reference Rf(kp) existe. 

Un nouveau vecteur d'erreur V^^ est alors nn6moris6 pour servir d'entr6e au traltement d6s points de la couronne 
suivante. 

20 Ainsi, on a defini un critere de proximite particulier selon lequel a partir du vecteur d'erreur entre les points d'un 

couple dans une couronne, on forme les couples dans la couronne suivante plus eloignee du centre, et. on mesure le 
nouveau vecteur d'erreur qui a grand! dans cette couronne suivante Et a partir du vecteur d'erreur mesure dans la 
couronne precedente, ainsi de proche en proche, la procedure est repetee jusqu'a ce que la taille des couronnes soil 
^gale ^ la taille de I'image. 

25 

ll-D/ CONSTRUCTION D'UNE GRILLE THEORIQUE PAR UN PROCEDE P'lTERATION FQNDE SUR L HYPOTHE- 
SE DES DISTORSIONS RADIALES (FIG.2, blocs 43 a 46> 

La methode selon ['invention consiste, h partir des caracteristiques de la premiere grille th6orique TG<^, h construire 
30 malntenant la meilleure grille theorique possible, c'est-a-dire a construire une grille theorique TG" aussi bien corrigee 
en distorsion que possible. 

Pour resoudre ce probleme la methode selon I'invention est iterative. Achaque iteration, cette methode cdmpren- 
dra les etapes de : 

35 * calcul d'un polynome F" de correction radiale, mis en oeuvre dans le bloc 44), 

* calcul de I'erreur E" de modelisation et memorisation de Toptimum, mis en oeuvre dans le bloc 45, 

* modification du centre optique OC^ mise en oeuvre dans le bloc 46, 

* estimation d'une nouvelle grille theorique TG", mis en oeuvre dans le bloc 43. 

40 Dans les etapes precedentes, on avait affecte d'un indice zero (0), place en exposant, toutes les notations relatives 

aux premieres estimations du centre GC^, du pas GMO et des couples P^ de points de la premiere grille theorique TG^ 
a construire. Chaque iteration sera reperee par un indice place en exposant, cet indice croissant a partir de zero a 
chaque nouvelle iteration. Par exemple, a IMteration n, le centre de grille est GC", le centre optique est DC", le pas de 
grille est GM", les couples sont P[J etc, 

45 En r6f6rence avec la FIG. 2, cette methode iterative forme une boucle. 

En supposant conn us les elements de construction d'une premiere grille theorique telle que TG^, d6finie par son 
centre de grille GC^, son pas GM°, son centre optique OC°, 

* on catcule d'abord, dans le bloc 44, un polynome F*^ de correction des distorsions radiales, 
so * on calcule ensuite dans le bloc 45.. une erreur E^ de modelisation, 

* on modifie alors, dans le bloc 46 fermant la boucle, la position du centre optique approximatif qui etait OC^ dans 
la premiere grille theorique et que Ton deplace en OC^ pour chercher a le rapprocher d'une position ideale ou 
rhypothese radiale est verifiee, c'est-a-dire oij les points de grille theorique sont alignes avec les points de refe- 
rence et avec le centre optique^ comme montre sur la FIG. 3D, 

55 * ayant modifie le centre optique maintenant note DC, on reestime dans le bloc 43 les elements de construction 
d'une seconde grille theorique TG^ definie par son centre de grille GC^ son pas de grille GM"" ; on construit cette 
seconde grille theorique TG'' au moyen de ces nouveaux element GC^ et GM^ puis toujours dans ce bloc 43, on 
realise la mise en correspondance des nouveaux points de grille theorique TRV^^^ et des points de r6f6rence de 
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rimage source SG^, 

* a ce moment la, on poursuit le procede d'iteration dans le bloc 44 par le calcul du polynome de correction des 
distorsions radiales, puis dans le bloc 45 par le calcul de I'erreur de modelisation et on reboucle I'iteration dans 
le bloc 46 par la modification du centre optique qui devierit OC^, etc. 

5 

On va se placer maintenant dans I'iteration n, en sortie du bloc 46 de la boucle, ou Ton vient de modifier la position 
du centre optique, qui devient OC". 

On porte alors cette connaissance du nouveau centre optique OC" dans le bloc 43 de la chaine d'iteration ou Ton 
dispose encore seulemen::des caracterlstiques de la grille theorique TG""^ k iteration n-1, et ou Ton va calculer des 
io nouvelles caracteristiques de grille theorique, un nouveau centre de grille GC" et un nouveau pas de grille GM" pour 
construire la nouvelle grille theorique a I'iteration n. 

Dans ce meme bloc 43, on effectuera ensulte la mise en correspondance des nouveaux points TR^^^^^ de grille 
theorique et des points de reference Rf^^j pour former des nouveaux couples PjJ. 

^5 II.D1/ Estimation precise du centre de grille GC". et du pas de grille GM " (FIG.2, bloc 43) 

En reference avec la FIG. 2, dans le bloc 43, a I'iteration n, on part de la connaissance du centre optique OC", 
nouvellement mis-^-jour dans le bloc 46 de la boucle d'it6ratlon. 
On part egalement des couples de points PjJ^ , qui sont not6es : 

20 

ou les Indices en exposant indiquant le numero de I'iteration n-1 ou ces valeurs ont ete determinees. 
On part aussi de la connaissance a ta derniere iteration n-1 : 

25 * du centre de grille GC"-! = (gcx"-"", gey"*"'), 

* du pas de grille theorique GM"-^ = (gmxf^-i, gmy"-i) 

On cherche a determiner maintenant k riteration n : 

30 * un centre de grille ameliore GC" = (gcx"^, gey") 

* et un pas de grille ameliore GM" = (gmx", gmy"), 

Ces parametres definissent une nouvelle mise-a-jour des points de grille TR^^^^. Toutes les caracteristiques de grille 

mises a jour sont calculees en pixels dans le repere de coordonnees fixe Ox, Qy. 
35 A cet effet, on utilise I'hypothese de base appelee hypothese radiale selon laquelle. si le centre optique OC" est 

correct, alors les distorsions qui sont radiales par nature, ont pour effet de deplacer les points de la grille de la mire 

reelle M, le long de rayons passant par le centre OC" et par les points de reference R^^^) comme montre sur la FIG. 3D. 
Cela indique que les points TR^^^^ de la grille theorique devraient se trouver pres de ces rayons. C'est-a-dire que 

la grille theorique la meilleure doit etre calculee de maniere a minimiser la distance geometrique entre les points de 
40 la grille theorique TR^J^^ a I'iteration n-1 et un rayon passant par le centre optique OC^ estime dans cette boucle 

d'iteration n et lesdits points de reference Rf^y^y 

En reference avec la FIG.7A, on a represents le centre optique OC" venant du bloc 46, le point de grille TR^jJ^ 

memorise a I'iteration precedente n-1 , le point de reference correspondant R^^,^) formant a I'iteration precedente n-1 

un couple P^'^ avec ledit point de grille TR^-J^. On a trace le rayon OC" - R^^j, et la distance geometrique AjJ""" entre 
45 le point de grille TR^'J^ et le rayon OC" - Rf^^^y 

Dans cette etape, k I'iteration n, on cherche k minimiser cette distance geometrique AjJ*^ , pour ameliorer I'aligne- 

ment du centre optique OC^ et des points des couples form6s des points de grille et des points de rSfSrence corres- 

pondants. 

Plus precisement, on definit une distance 



so 



55 



Ak=[TR9;K,.(OC".R„,))] 

qui est la distance geometrique memorisee a I'iteration n-1 , entre un point TR^'j^ de ta grille theorique TG"-i et le rayon 
passant par le centre OC" et le point de reference correspondant R^^^^. 
On peut ecrire que, dans le repere Qx, Oy fixe : 

__n-1 r/ n-1 - . , n-1 n-1 

TRg^k)=t(9cx +gmx x (gey +gmy x ^ly^^^)] 

Dans cette formule les points de grille TR^'J^ ne dependent que des composantes gmx'^'^ gmy"-^ du pas de grille 
theorique et des composantes gcx"'\ gcy""^ du centre de grille theorique. 



14 

BNSOOCID: <EP ^07201 25AlJ_> 



EP 0 720 125 A1 

Dans cette formule §galement, les param^tres et jiy^,^) d6finissent respect ivement le nombre de pas de grille 
th6orjque respeclivement gmx"'"" et gmyV"'' dont il faut se deplacer parallelement aux axes Qx, Qy, pourtrouver le point 
de grille TRgJ,[j a partir du centre de grille theorique de coordonnees gcxf^'^ gey"'"*. 

Par exemple. pour realiser la FIG.7A, on a suppose que le point de grille etait a X^^,^) = 4 pas de grille en horizontal 
s elk My^,^) = 2 pas de grille en vertical a partir de grille GC"-^ lui-meme rep^re par ses composantes gcx"-"". gcy"-^ en 
pixels dans le rep^re fixe Qx, Oy. 

Ensuite, le centre GC" et le pas de grille GM" sont calcules a cette iteration n dans le bloc 43 en effectuant la 
minimisation d'un critere note O**, appele ci-apres critere radial, qui exprime I'hypothese radiale enoncee plus haut, 
seton laquelle la meilleure grille theorique possible correspond d des distances A[] plus petites k Titeration n qu'^ 
10 iteration n-1. 

Ce critere radial s'dcrit : 

$"(gcx\ gcy^, gmx", gmy") = 2 (Aj)^ 

k 

Pour que ce critere <E»" ait une solution non ambigue, il convient d'exclure la solution ou le pas de grille est nul et 
ou le centre de grille estconfondu avec le centre optique, car dans tous ces cas le critere est nul, et aucune solution 
constructive ne peut etre trouvee. Par contre. si le pas de grille est contraint d'avoir une grandeur raisonnable, et si le 
centre de grille est different du centre optique, alors le probleme est bien pose. 

Toute m6thode de minimisation, menant k une solution non singulifere. est bonne et peut §tre utilis6e. On propose 
selon r invention une methode incluant les stapes suivantes : 

* depart avec une valeur raisonnable du pas de grille GM"""", en posant dans le critdre radial 4>", les composantes 
gmx"-i et gmyv-i du pas de grille fixes, 

* minimisation du critdre radial O", par rapport au centre de grille GC"-"* : les composantes gmx"^""", gmy^""" du pas 
etant fixes, on derive la formule du critere radial ^ par rapport aux deux coordonnees variables gcx"""", gcy""i du 
centre de grille GC^'V Ce calcul de derivees fournit deux formules expllcites qui donnent des valeurs de coordon- 
nees gcx" et gey" pour mettre a jour le centre de grille GC" a cette iteration n, 

* minimisation du critere radial <I>", par rapport au pas de grille GM*^-! : on fixe maintenant les coordonnees gcx" et 
gey" du centre de grille qui viennent d'etre calculee, et on derive la formule du critere radial par rapport aux 
composantes du pas de grille variables gmx^^^i, gmy""^ Ce calcul de derivees fournit a nouveau deux formules 
explicites qui donnent des valeurs de composantes gmx" et gmy^ pour mettre a jour le pas de grille GM" a cette 
iteration n, 

* repetition de ces operation autant de fois que necessaire jusqu'a I'obtention d'une precision requise. 

55 li-D2/ Estimation a I'lteratlon n des couples des points de reference-points de grille theorique correspondants 
(FIG.2, bloc 43) 

Cette operation est faite en pratiquant exactement la meme procedure que pour la premiere grille theorique k 
I'exception que Ton utilise les valeurs du centre et du pas GC" et GM", respeclivement mises k jour k I'it^ration n lors 
40 ' de I'etape precedente. 

A I'iteration n, dans le bloc 43 de la FIG. 2, on dispose maintenant dans le repere fixe Qx, Qy, de ; 

* une mise k jour des coordonnees du centre de grille GC" = (gcx", gey") 

* une mise a jour des composantes du pas de grille GM" = (gmx", gmy") 
45 * la position du centre optique OC" provenant du bloc 46. 

Ces caract6ristiques permettent de construire une nouvelle grille th6orique TG" mise k jour, k mailles rectangu- 
laires regulierement espacees, comme sur la FIG.3E. 

On realise alors la procedure de mise en correspondance des points de r6f§rence Rf^|^^ de la grille source SG*^ et 
50 des points '^Gg(k) de la nouvelle grille theorique TG*^ pour former des nouveaux couples P^. En effet, du fait du chan- 
gement du centre et du pas de la nouvelle grille theorique, le point de reference qui avait ete trouve k I'iteration n-1 
pour former un couple avec un point de grille donne peut maintenant ne plus etre le meilieur candidat. Une nouvelle 
mise en correspondance doit done etre effectuee. 

On rdalise alors cette procedure de nouvelle mise en correspondance des points de grille de la grille theorique 
55 TG" et des points de reference de I'image source en considerant des couronnes rectangulaires croissant autour du 
centre de grille GC" localise par ses nouvelles coordonnees (gcx", gey"), et en formant les couples P[J par progression 
pas-a-pas, couronne apres couronne, en appliquant le critdre de proximity propose dans I'expose de la construction 
de la premiere grille theorique TG^. 
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ll-E/ CALCUL D'UN POLYNOME DE CORRECTION RADI ALE (FIG.2, bloc 44) 

En partant des couples P[J determines a I'etape precedente de la chaine qui constitue I'iteration n. on calcule des 
couples de rayons : 

* un premier rayon X[J ayant pour valeur la longueur du segment qui lie le centre optique OC" et point de grille TG" 
d'un de ces couples PjJ ; 

* un second rayon YJJ ayant pour valeur la longueur du segment qui lie ie centre optique OC" et ie point de reference 
Rf(k) de ce couple PjJ . 

En reference avec la FIG.7B, ces rayons s'expriment par : 

<=(OC".TR^,,,) 

is ■ y;:=(oc".r„^,) 

Ces couples de rayons sent calcules pour chaque couple PjJ de la grille th§orique TG" nouvellement mise a jour. 
On calcule ensuite une fonction polynomial F" qui relie au mieux les rayons Y]J aux rayons XP de chaque couple 

20 En reference avec la FIG.8A, on a represente sur un graphe des rayons XjJ^ , XjJ^, en abscisse X, et des rayons 

YjJ^ , Yjjgen ordonnees Y, formant des couples PjJ^ respectivement, le centre optique OC" etant place a I'origine 
des axes de coordonnees graduees en pixels. Ces couples , PjJ^ forment, sur ce graphe, des points symbolises 
par des petits ronds. qui se repartissent de part et d'autre de la courbe representative de la fonction polynomiale F" 
recherchee. Sur la FIG.SAon s'est limite a la representation de deux couples de rayons pour la clarte du dessin. 

25 En reference avec la FIG. SB, on a represente sur un graphe dont les axes de coordonnees X et Y sent gradu6s 

en nombre de pixels, six couples OP^ , P^. P^. Pg^, pg^, P^g formes en relation avec la premiere grille theorique 
SGO. Sur ce graphe, ie centre optique 00° est porte a I'origine des axes de coordonnees, les rayons X^ , Y^ sont portes 
sur I'axe des abscisses et des ordonnees respectivement. A cette iteration numero zero, les points representant les 
couples P° forment un nuage de points autourde la courbe representative de la fonction F^ trouvee pour relier les 

30 rayons XJJ, Y[J 

En reference avec la FIG.8C, on a represente sur un graphe dont les axes de coordonnees X et Y sont gradues 
en nombre de pixels six autres couples P;;; . P;^ , PJ^g . PJ^^ , P|^ , P^^^ formes en relation avec une grille theorique SG*". 
k une iteration numerotees m ou la correction des distorsions est devenue pratiquement parfaite. Dans ce cas I'en- 
semble des couples constitue une courbe lissee correspondant a la fonction polynomiale F"^. 
35 En reference avec ces FIG.8A a 8C, on a trouve que la fonction F"^ la mieux appropriee a relier les couples de 

rayons X[J , Y[J a iteration n. etait un polynome de degre 5 ou 6 ; un polyndme de degre 3 par exemple ne suffit pas a 
resoudre ce probleme. II est apparu que la correction des distorsions radiales necessitait le calcul d'un polynome 
d'ordre superieur a 3. 

En reference avec la FIG.8A, a I'iteration n, on peut mesurer pour chaque rayon XjJ^ , XjJ^... Terreur respectivement 
40 E[J^ . EjJ^-- qui existe entre les rayons Y|J^ , YJJg... et la fonction polynomiale F" qui relie au mieux les points de couples 
' ^k2 - • ''ensemble des couples PjJ de la grille theorique TG", a I'iteration n, cette erreur appelee ci-apres erreur 
de modelisation E" s'ecrit : 

E" = 2 {[YiT - F^CXiT)]^} 

45 k 

11 est important de noter que le polyn6me F", cl-apres appele polyndme de correction radlale, constitue ^ 
riteratlon n, une formula de transformation qui permet de passer : 

* d'un point TR^^,^^ de la grille theorique corrigee, 

so * ^ un point de reference P{^y^ de la grille source distordue. 

Ce polynome F" peut etre decompose en deux composantes rectangulaires F" et F". 

X y 



55 



ll-F/ CALCUL DE L'ERREUR DE MODELISATION rFiG.2, bloc 45> 

En reference avec la FIG.2, dans la chaine qui constitue Titeration n, a I'etape symbolisee par le bloc 45, on calcule 
le meilleur polynome de correction radiale F" pour minimiser I'erreur de modelisation E", par une operation de regres- 
sion lineaire classique. 
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Une erreur de modelisatioo E" est calcul6e k chacune des iterations de 0 k, par exemple. 10. Et ces erreurs de 
modelisation E^, E^ .. E""^ sont memorisees. Parmi les differentes erreurs de modelisation calculees au cours des 10 
cycles d'lt6ration pris comme exemple, I'une d'elle. a une iteration nunnerotee Nopt. se revele inferieure aux autres 
erreurs calcul6es aux iterations ant^rieures et everitueilem^nt aiux iterations ulterieures. 
5 Cette erreur minimale E^opt est nnemorisee en nnenrje temps que le polynome de correction radiale pNopt, ou bien 

ses composantes rectangulalres f^^^\ f^opt 

X y 

l-G/ MODIFICATION DU CENTRE OPTIQUE DE LA GRILLE THEORIQUE (FIG.2, bloc 46) 

10 L'erreur de modelisation E" calcul6e pr6c6demment k I'it^ration n, dans la methode selon {'invention bas6e sur 

rhypothdse radiale, depend fortement de la position OC" trouvee pour le centre optique, qui devient ideal lorsqu'il est 

au centre des distorsions. 

Dans la chame d'iteration representee sur la FIG. 2, on trouve une boucle concretisee par le bloc 46, qui regoit en 
entree les donn^es du bloc 45 de calcul d'energie d'erreur et qui fournit en sortie, au bloc 43, un centre optique modifie 
IS 0C"+^ pour calculer une nouvelle grille theorique TG"+'' dans une nouvelle iteration n+1 . 

On rappelle que selon I'invention, la minimisation de I'energie d'erreur E" fonne la solution du probleme de cor- 
rection des distorsions fonde sur I'hypothese radiale. 

A cat elfet. ^ chaque iteration, de 0 a 10 par exemple, la position du centre optique OG", repere dans la grille 
th6orique TG", par ses coordonn6es dans ie rep§re fixe Qx, Qy 

20 

OC" = (ocx", ocy") 

est modifiee pour cooperer a minimiser l'erreur de modelisation E". La meilleure position du centre optique OC" pour 
obtenir cette minimisation peut etre recherchee de plusieurs manieres. 

Dans une premiere methode, une zone de recherche large est definie autour du centre optique de depart que Ton 
2S avait note OQO a Titeration zero. Au cours de chaque iteration, on teste un des points de cette zone de recherche large, 
pour constituer un centre optique. On realise autant d'iterations que de points a tester dans la zone de recherche large 
(FULL SEARCH). On retient finalement comme centre optique celui qui correspond a l'erreur de modelisation la plus 
faibie. 

Une seconde methode, moins couteuse en temps de calcul, estfondee sur le calcul de gradients. A deux iterations 
30 successives, une premiere localisation, puis une seconde localisation sont essayees. dans une direction donnee a 
partir du premier centre optique OC^. Selon que ces premiers essais condulsent a une diminution ou a une augmen- 
tation de l'erreur de modelisation E", la meilleure position du centre optique continue a etre recherchee dans la.meme 
direction et le meme sens, ou bien dans la meme direction et en sens oppose. Puis d'autres directions sont essayees 
. faisant un angle donne avec la premiere direction et finalement on retient Tlt^ration ou, dans une direction et un sens 
55 donne, on a trouve le meilleur centre optique pour minimiser l'erreur de modelisation. 

On memorise le numero NIast de I'iteration a laquelle on a trouve la position de ce centre optique la meilleure pour 
minimiser l'erreur de modelisation. 

il-H/ SORTIE DE LA CHAiNE D'ITERATION 

40 

Dans ces conditions les meilleures conditions de mise a jour en sortie du bloc 45 represents sur la FIG.2 corres- 
pondent a une iteration ou Ton fait cooperer les conditions trouvees a Nopt et les conditions trouvSes a NIast. Cette 
iteration devient la derniere iteration : elle est numerotee Nopt (NIast). 

La meilleure fonction polynomlale correspondante, pour minimiser l'erreur de modelisation, est alors notee : 

Q ^ pNopt(Nlast) 

et ses composantes rectangulalres sont deux polyndmes notes : 

Nopt (NIast) 

so Gx = ^y) 

Nopt (NIast) 

ss Ces composantes sont des polyndmes ^ deux coordonnees car chacun doit etre applique k un point couple P[J en 
coordonnees cartesiennes. 

En reference avec la FIG.7C et avec la FIG. 3D, il apparait qu'apres les operations d'iterations qui ont men6 k la 
determination des conditions correspondant a la derniere iteration NOpt (NIast) le centre optique a nnaintenant une 
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position : 

OC^last 

et les points de grille TR^|^* et de reference Rf^j^j formant les couples p^^ast qq^^ align6s au mieux sur un mdme rayon. 
En reference avec la FIG.8D, les points de couples p^^^^ sont sur la courbe lissee repr^sentant la fonction poly- 
s nomiale pNiastCNopt) 

line grille theorique ideale correspondant a ces conditions est montree sur les FIG.3F et 3G. TG'^"' et TG"^ 
Une fois que la fonction polynomiale G, ou bien les composantes Gx et Gy ont ete calculees, TOUTES LES ETAPES 
PRECEDENTES SONT TERMINEES DEFINITIVEMENT 

La connaissance des fonctions G, ou Gy, Gy est suffisante pour assurer la correction des distorsions optiques de 
10 toute image enregistr^e par la camera 2 et son syst^me de prise de vue 1 (FIG.1). 

III. CORRECTION DES DISTORSIONS OPTIQUES DANS UNE IMAGE (FIG.1, sous-ensemble 47 : et FIG.2. sous- 
ensemble 47) 

Comme on I'a dit precedennment, les composantes Gx et Gy sont des formules de transformation qui font passer 
de I'espace rectifie note image cible 71, vers I'espace distordu note image source SI d'origine. 

Dans rimage cible k construire Tl, les coordonnees de chaque pixel Xj, et y-j-i dans le repere fixe iix, r2y, sont 
memorisees. L'image Tl est alors balayee pixel par pixel, par exemple par balayage video usual 

La valeur de Tintensit^ ou niveau de lumi^re monochromatique, ou niveau de gris qui doit etre attribu6 au pixel 
ecu rant Xj,, yj, de ('image cibie, est cherchee a une adresse de I'image source X31, ygj calcul6e k I'aide des polynomes 
de transformation Gx, Gy stockes dans le bloc 47 : 

Xs» = Gx(XT|.yT!) 

Le point trouve a I'adresse Xsi, ysi dans I'image source SI n'etant en general pas un pixel, le niveau de lumiere 
monochromatique ou niveau de gris est calcule pour ce point par interpolation a partir de I'intensite de ses voisins 
pixels dans cette image source. Puis la valeur d'intensite calculee par interpolation est attribute au pixel courant k 
I'adresse x-pi, yji dans I'image cible et est memorisee dans la seconde memoire d'image. 



Revendlcations 

1. Dispositif de formation d'image comprenant ; 
un systeme de prise de vue (1), 

un systeme de saisie (2) des donnees d'une image source (SI) distordue par le systeme de prise de vues, 
un systeme de numerisation (3) de ladite image source (SI) incluant une premiere memoire d'image pour le 
stockage des donnees d'intensite de chaque pixel repere par une adresse dans une matrice bidimensionnelle, 
un systeme de traitement d'image numerique (4,47) de construction d'une image cible (Tl) corrigee en distor- 
sions, correspondant a I'image source (SI), 

dans lequel, le systeme de traitement d'image inclut : 

un premier sous-ensemble (4) de determination prealable a partir des donnees d'image d'une image source 
TEST distordue (SG-2), d'un centre optique (OC") et d'une fonction polynomiale (F Nopt/NLast) cle correction 
des distorsions radiales autour de ce centre optique. 

un second sous-ensemble (47) incluant une memoire de stockage du centre optique et de la fonction polyno- 
miale determines par le premier sous-ensemble, un bloc de calcul qui applique la fonction polynomiale de 
correction a chaque adresse de pixel de I'image cible (Tl) pour fournir I'adresse d'un point dans rimage source 
distordue (SI) ou trouver une donnee d'intensite a appliquer a I'adresse de depart dans I'image cible, et une 
seconde mdmoire d'image pour le stockage des donnees de I'image cible (Tl) reconstruite. 

2. Dispositif de formation d'images seion la revendication 1 , dans lequel: 
55 

le premier sous-ensemble (4) inclut des blocs de traitement de {'image source TEST, ladite image source 
TEST etant en forme de grille appelee grille source (SG^) pour construire une image cible TEST appelee grille 
theorique (TG"). lesdits blocs comprenant : 
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un bloc d'extraction (4.1 ) de points de rSf^rence (R|^) aux intersections des barreaux de la grille source (SG^), 
un bloc d'estimation (42) d'une premiere grille theorique (TG^) approchee, calculant I'adresse d'un centre 
(GC^) et ia grandeur d'un pas (GM^) de ladite premiere grille theorique (TG^), 

un bloc d'estimation (43) d'une grille theorique (fG") dorrigee en distorsion, calculant par iteration I'adresse 
d'un centre (GC") et la grandeur d'un pas (GM") de ladite grille theorique corrig^e (TG"), 
un bloc de calcul (44) d'un polyndme de correction des distorsions radiales fournissant une toi de transforma- 
tion pour passer d'un point de la grille theorique corrigee a I'iteration n (TG") k un point de r6f6rence (R^) de 
la grille source distordue, loi de transformation a I'iteration n, 

un bloc de calcul (45) d'une erreur de modelisation (E") due k la loi de transformation k I'iteration n, 

un bloc de calcul (46) d'une adresse modifiee du centre optique (Oc") qui minimise I'erreur de modelisation. 

Precede pour effectuer des corrections des distorsions optiques geometriques produites dans une image par un 
systeme de prise de vue, ce precede comprenant des etapes preliminaires de : 

a) saisie des donnees et numerisation d'une image source TEST (SG^), incluant le stockage des donn^es 
d'intensite de chaque pixel repere par une adresse dans une matrice bidimensionnelle, 

b) estimation d'un centre optique (OC") localise au mieux au centre des distorsions optiques de I'image source 
TEST (SG^) et report de ce centre optique dans une image numerique a construire (TG"), referencee image 
TEST cible, repr6sentant I'image source TEST corrig6e des distorsions optiques, 

c) estimation d'une fonction polynomiale (F") pour faire correspondre a I'adresse d'un pixel (TR[J^) dans I'image 
TEST cible corrigee, I'adresse d'un point correspondant {H^), note point de reference, dans I'image source 
TEST distordue, en se fondant sur I'hypothese de base que les distorsions optiques geometriques dans I'image 
source TEST sont radiales autour du centre de distorsion, de sorte que lesdits points correspondents (TRjJ^, 
R|^) sont alignes au mieux avec le centre optique estime (OC"). 

Procede selon la revendication 3, comprenant en outre dans les etapes preliminaires des etapes iteratives de : 

d) .estimation de la meilleure fonction polynomiale (F^opt) capable de minimiser, a une iteration d'ordre Nopt, 
I'erreur de modelisation (E") realisee sur la localisation des pixels de I'image TEST cible (TG") construite k 
une Iteration precedente (n) selon I'hypothese des distorsions radiales autour d'un centre optique (OC") estime 
a ladite Iteration precedente. 

e) estimation, a une iteration d'ordre NIast, d'un centre optique modifie optimise (OC^'^®^) et reestimation dans 
ces conditions d'une nouvelle fonction polynomiale (pNopKNiast)) capable de minimiser. davantage I'erreur de 
modelisation realisee sur la localisation des pixels de I'image cible reconstruite selon I'hypothese des distor- 
sions radiales autour dudit meilleur centre optique (OC^'^^) estime. 

Procede selon la revendication 3, dans lequel les etapes preliminaires sontfaites une fois pourtoutes, ce proc§d6 
comprenant a Tissue des etapes preliminaires, des etapes de : 

aj saisie des donnees d'une image source (SI) quelconque distordue par le dispositif de prise de vue, et 
numerisation et stockage des donnees d'intensite de chaque pixel, 

f) correction proprement dite des distorsions optiques, en appliquant a une image cible numerique (Tl) a cons- 
truire, ia fonction polynomiale (F") estimee, de maniere a fournir a partir de chaque adresse de cette image 
cible (Tl). une adresse correspondante dans ['image numerique source (SI) distordue ou trouver une donn6e 
d'intensite k appliquer k ladite adresse dans I'image cible k construire. 

Procede selon la revendication 4, dans lequel les etapes preliminaires sontfaites une fois pour toutes, ce precede 
comprenant a Tissue des etapes preliminaires, des etapes de : 

aj saisie des donnees d'une image source (SI) quelconque distordue par le dispositif de prise de vue, nume- 
risation et stockage des donnees d'intensite de chaque pixel ; 

f) correction proprement dite des distorsions optiques, en appliquant, a une image cible numerique (Tl) k 
construire, la fonction polynomiale (Ff^opKNiast)) estimee la meilleure, de maniere a fournir a partir de chaque 
adresse dans cette image cible (Tl), une adresse correspondante dans I'image numerique source (SI) distor- 
due ou trouver une fonction d'intensite k appliquer k ladite adresse dans I'image cible a construire. 

Precede selon Tune des revendications 3^6, selon lequel : 
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dans r^tape a) preliminaire de saisie de donnees et numerisation, Timage source TEST est I'image numSrisee 
d'une grille a mailles rectangulaires dont ies barreaux sont parall^les aux lignes et colonnes de pixels, aux 
erreurs de distorsions optiques pres, referencee grille source, 

dans Tetape b) I'estimation d'un centre optique (OC") est effectuee de maniere iterative a partir d'une iteration 
5 zero comprenant la construction d'une premiere image TEST cible corrigee en distorsion, notee premiere grille 

cible (TGO) au moyen des sous-etapes de : 

extraction d'un point par zone d' Intersection des barreaux de la grille source, ces points extralts etant 
not6s points de reference (R^) de la grille source (SG^), 
10 estimation d'un point de reference (R^) le plus proche du centre de distorsions de la grille source, 

report de ce point de reference (R^) dans la premiere grille cible (TG^) pour constituer le premier centre 
de grille (GCO), 

estimation d'un premier pas de grille (GM^) pour iadlte grille cible (TG^): 

estimation d'un premier centre optique approximatif (OC^) pour ladite premiere grille cible (TG^), situd en 
IS coincidence avec ledit point de reference (R^) reporte et en coincidence avec le centre de grille (GG°). 

8. Proc6d6 selon la revendication 7. dans lequel a t'etape d) Testimation iterative de la meilleure fonctlon polynomiale 
(pNopt) comprend I'estimation a I'iteration d'ordre zero d'une premiere fonctlon polynomiale (F^) comprenant Ies 
sous-6tapes de : 

20 

construction des points (TR^) de ladite premiere grille cible a partir du centre de grille (GC°) et du pas de 
grille (GM^) pour correspondre aux points de reference (Ri^) d'intersection des barreaux de grille source (SG^) 
mise en correspondence des points de grille (TR^) de la premiere grille cible (TG^) et des points de reference 
(R|^) de la grille source (SG^) en localisant le centre de grille (GC^) en coincidence avec le point de reference 
25 (rO) le plus proche du centre de distorsion, 

formation des couples constltues des points de grille (TR^) et des points de reference (R,J correspondants, 
de proche en proche a partir du centre vers Ies bords de la grille cible (TG^), 

estimation de couples de rayons (X^, Y^) formes de segments joignant le centre optique (OC^) de la premiere 
grille cible (TG^) et chaque points des couples points de grille cible-point de reference (TR^, R,^), 
30 caicul de la premiere fonction polynomiale (F^) comme celle qui relie le mieux Ies couples de rayons (X^, Y^) 

pour respecter I'hypothese selon laquelle Ies distorsions sont radlales autour du centre optique connu (OC^). 

9. Precede selon la revendication 8, dans lequel a I'etape d) I'estimation Iterative de la meilleure fonction polynomiale 
(pNopt) comprend I'estimation a une Iteration d'ordre n > 0, d'une fonctlon polynomiale (F") comprenant Ies sous- 
es etapes de : 

construction d'une grille cible (TG") a t'ordre n, ayant un centre optique (OC") et des points de grille (TR[J,) 

definis par un centre de grille (GC") et un pas de grille (GM"), 

mise en correspondance des points de grille (TRfJ^) de la grille cible (TG") et des points de reference (F^) de 
40 la grille source (SG") en localisant le centre de grille cible (GC") par ses coordonnees determjnees a I'iteration 

precedente, 

formation de couples (PJJ) de points de grille (TR[J^)-polnt de reference (l\) de proche en proche a partir du 
centre vers Ies bords de I'image cible (TG"), 

estimation des couples de rayons (X[^, Y[J) formes des segments joignant le centre optique connu (OC") et 
45 chaque point des couples (TRJJ^, R^), 

caicul de la fonctlon polynomiale (F") comme celle qui relie le mieux Ies couples de rayons (XjJ, YjJ) pour 
respecter Thypoth^se selon laquelle Ies distorsions g6om6triques sont radlales autour du centre optique (OC"). 

10. Procede selon la revendication 9, dans lequel la construction d'une grille cible (TG") a I'iteration courante (n) 
. ^ comprend Ies sous-stapes de : 

localisation d'un centre optique (OC") nouvellement mis-^-jour, 

localisation des points de grille (TR[J^'^ ) et points de reference (R,^) correspondants constituant des couples 
(PjJ"^ ) formes a I'iteration precedente (n-1), 
55 estimation des distances geometriques (A[J*^ ) abscisses des points de grille (TRU[^ ) sur Ies rayons passant 

par le centre optique nouvellement mis a jour (OC") et par Ies points de reference 
correspondants (R|^) des couples (PJ'^ ) formes a cette ItSratbn precedente (n-1), 

estimation d'un critere reference critere radial qui exprime I'hypothese radiale que la meilleure grille cible 
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corrigee correspond ^ la minimisation des distances geometriques (AjJ*"* ), 

minimisation du critere radial exprime en fonction des coordonnees du centre et des composantes du pas 
de la grille a I'ileration (n-l), fournissant Jes cpordonnees du centre et les composantes du pas de la grille 
cible k I'iteration n. ^ - 

5 

11. Precede selon la revendication 10, dans lequel I'etape e) d'estimation du meilleur centre optique (OC^'ast) est 
realisee dans une boucle d'iteration dans laquelle a chaque iteration courante (n) la localisation du centre optique 
(OC") est modifiee, puis reintroduite dans la construction de la grille cible courante (TG") pour la determination 
d'une nouvelle fonction polynomiale (F") et I'estimation de I'erreur de modelisation (E") correspondante, la meilleu- 
10 re localisation du centre (OC'^'ast) ©tant celle qui resulte en I'erreur de modelisation la plus faible correspondant d 

la meilleure fonction polynomiale (FNopt(Niast)) 
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